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ABSTRACT: The use of Korean native chicken is increasing, and the discovery of new genetic resources is very important from both 
economic and genetic conservation points of view. In this study, mtDNA D-loop sequences from 272 privately-owned Korean native 
chickens from a Hyunin farm were investigated. Seventeen nucleotide substitutions were identified from the sequence analysis and they 
were classified as 6 haplotypes. Previously investigated haplotypes in five Korean native chicken populations have been compared with 
the Hyunin chicken population. The results indicated that two haplotypes, HIO and H15, in the Hyunin chicken population were not 
previously identified in other Korean native chicken populations, representing 33.09% (90/272) and 1.1% (3/272) of the Hyunin 
population, respectively. On the other hand, four other haplotypes were identical to those of a previous study of Korean native chicken 
populations. This result is indicative of conservation strategies of Hyunin chicken populations for expanding the genetic diversity in the 
Korean native chicken population. (Key Words: D-loop Variation, Haplotype, Korean Native Chicken, mtDNA) 



INTRODUCTION 

Increased research and genetic profiling of native 
species is an important component of modem perspectives 
of genetic diversity conservation. There is also an 
increasing demand for native chicken meat currently, even 
though 90% of chicken consumed in Korea is from 
imported breeds (MIFAFF, 2009). Five lines of Korean 
native chicken populations were developed in the 1990's, 
and this breeding was mainly controlled by a national 
institution called National Institute of Animal Science 
(NIAS) in Korea. However, there are still native chicken 
breeds available in the private farms, and one of them, the 
Hyunin population, is used in this study for comparison 
with the primary Korean native chicken populations for 
which we have previous genetic information. 

Maternally-inherited mitochondrial DNA (mtDNA) 
sequences have been successfully used for genetic diversity 
studies of avian and domestic animals (Komiyama et al., 
2003, 2004; Liu et al., 2006; Odahara et al., 2006; Sasazaki 
et al., 2006; Lei et al., 2007; Wang et al., 2007; Li et al., 
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2008). In the mitochondrial genome, coding genes are 
predominantly used for biodiversity studies of many 
different species, whereas, the non-coding control region is 
used for the inter- specific population studies (Baker and 
Marshall, 1997; Moore and Defilippis, 1997). MtDNA has 
the advantages of the absence of recombination and rapid 
nucleotide substitutions compared with genomic DNA 
sequences (Aquadro and Greenberg, 1983; Lansman et al., 
1983; Cann et al., 1984). Chicken D-loop sequences were 
previously determined to investigate the origin of breeds or 
species, and various haplotypes were reported (Komiyama 
et al., 2003, 2004; Liu et al., 2006; Oka et al., 2007). 
Previously, three SNPs of D-loop region were found to have 
significant discrimination power for two commercial native 
chicken populations in Korea (Hoque et al., 2011). Also, the 
D-loop hypervariable region of mtDNA was used to 
determine the relationships between Korean native chickens 
and other chicken breeds (Hoque et al., 2009). 

In this study, the chicken mtDNA D-loop region was 
used for the demarcation of haplotypes in the privately- 
owned Hyunin native chicken population, and was 
compared with those of other Korean native chicken breeds, 
which were previously developed in NIAS. This Hyunin 
population has distinct feather colors and they were 
classified as eleven colors. Some of these (red, yellow, 
black, gray and white) are the same colors identified in 
Korean native chickens at NIAS. These results provide a 
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basis for an appropriate conservation breeding program for 
the Korean native chicken. 

MATERIALS AND METHODS 

Sampling and DNA extraction 

A total of 272 native chicken samples were collected 
from a single flock in a Hyunin chicken farm. Genomic 
DNAs were extracted from blood samples using a manual 
extraction method (Miller et al., 1998) applied to a large 
volume, which was stored at -20°C until use. For the 
comparison of haplotypes, previous D-loop sequence data 
from two commercial chicken populations Hanhyup (40 
samples) and Yelhm (37 samples) and three Korean native 
chicken breeds (88 samples) were included in the analysis 
(Hoque et al., 2011). 

PCR amplification and sequencing 

PGR was used to amplify a 591 bp fragment in the 
D-loop hypervariable region in mtDNA by the following 
primer pair (Forward: 5'- AGGAGTAGGGGTTGAAAAGG 
-3' and Reverse: 5'-ATGTGGGTGAGGGAGGAAGGAG - 
3'). The PGR reactions included approximately 100 ng of 
genomic DNA, 2.5 |il of lOx buffer (Tris-HGl (pH 9.0), 
PGR enhancers, (NH4)2S04, 20 mM MgGlj), 2.0 )il of 10 
mM dNTPs mixture (2.5 mM each of dATP, dGTP, dGTP 
and dTTP), 1 |il of 10 pM of each primer and 1 U HS Prime 
Taq (GeNet Bio, Korea) in a 25 ^1 reaction volume. PGR 
was performed in a My-Genie96 Thermal Block (Bioneer) 
with an initial denaturation step at 94°G for 10 min, 
followed by 35 cycles of 30 s at 94°G, 30 s at 61°G, 40 s at 
72°G and a final extension step at 72°G for 10 min. The 
PGR products for mtDNA were electrophoresed on 1.5% 
agarose gels with ethidium bromide, and DNA bands were 
visualized under ultraviolet light. Purification of PGR 
products was performed using an Accuprep®PGR 
purification kit (Bioneer) according to the manufacturer's 
instructions. Purified PGR products were also confirmed 
using agarose gels for sequencing. All the purified PGR 
products were sequenced by Genotech (www.genotech.co.kr) 
commercial company. 



Analysis 

The chicken mtDNA D-loop nucleotide sequence data 
were aligned using the GlustalW program (Thompson et al., 
1994). Nucleotide replacement export data from mtDNA 
were carried out in haplotype sequences by using MEGA5 
(Tamura et al., 2011). The number of haplotypes, nucleotide 
variable site, haplotype diversity and nucleotide diversity 
(Nei, 1982) were calculated using DnaSP V.5.10 (Rozaset 
al., 2003). A rectangular Kimura 2-parameter model 
neighbor-joining (NJ) phylogenetic tree with 1000 bootstrap 
replications was constructed using MEGA5. Also, an 
outUne for median-joining network profiles with positional 
SNPs for the differentiation of haplotypes by using 
NETWORK4610 software (Bandelt et al., 1999). 

RESULTS AND DISCUSSION 

Genetic diversity and haplotype analysis 

Mitochondrial D-loop sequence variations were used to 
calculate the number of variable sites (S), number of 
haplotypes (h), haplotype diversity (Hd) and nucleotide 
diversity (Pi) among the chicken populations (Table 1). 
Analysis of 272 mtDNA D-loop sequence data from the 
Hyunin population revealed six haplotypes, which were 
comprised of seventeen nucleotide substitutions. All of the 
six mtDNA haplotypes for the Hyunin chicken population 
that were discovered are shown in Table 2. Based on the 
statistical analysis, the identified Hd and Pi were 0.768 and 
0.0083, respectively, in Hyunin chickens. Based on the 
previous analysis of 40 sequences from Hanhyup 
commercial chicken population, 15 variable sites were 
identified which have 5 haplotypes with 0.778 of Hd and 
0.0068 of Pi. Interestingly, 37 mtDNA D-loop sequences 
were conserved in the Yellim commercial chicken 
population (Hoque et al., 2011). Moreover, three Korean 
native chicken breeds have been compared with the three 
Korean native chicken populations, dominated by black, red, 
and yellow, respectively. The black Korean native chicken 
(KB) breed was observed to be highly variable in terms of S, 
h, Hd and Pi for 22, 9, 0.823 and 0.0113, respectively Red 
Korean native chicken (KR) breeds also had a high 
nucleotide substitution of 20, and 6 haplotypes which used 



Table 1. The calculated number of variable sites (S), number of haplotype (h), haplotype diversity (Hd) and nucleotide diversity (Pi) 
using mitochondrial D-loop sequence polymorphisms for the Korean native chicken populations 



Breed/population Sample no. S h Hd Pi 



Hyunin (Hn) 


272 


17 


6 


0.768 


0.0083 


Hanhyup (Hf) 


40 


15 


5 


0.778 


0.0068 


Yellim (Yf) 


37 


0 


1 


0 


0 


Black Korean native chicken (KB) 


28 


22 


9 


0.823 


0.0113 


Red Korean native chicken (KR) 


30 


20 


6 


0.662 


0.0093 


Yellow Korean native chicken (KY) 


30 


10 


4 


0.591 


0.0077 



S = Number of variable site.s, h = Number of haplotype, Hd = Haplotj'pe diversity, Pi = Nucleotide diversity. 
Except for the sequences from Hyunin, the sequence data were obtained from Hoque et al. (2011). 
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Table 2. The identified haplotypes witli tlie mitochondrial D-loop sequence polymorphisms in Hyunin chicken 



Nucleotide substitution^ 





1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


2 


3 


3 


3 


3 


3 


4 




Haplotype 


6 


9 


1 


1 


1 


2 


4 


4 


5 


6 


9 


1 


1 


1 


6 


7 


4 


Sample no. 




7 


9 


0 


2 


7 


5 


3 


6 


6 


1 


1 


0 


3 


5 


7 


0 


6 




Hnhapl 


C 


T 


c 


G 


T 


T 


T 


C 


T 


C 


G 


c 


T 


C 


C 


c 


C 


29 


Hnhap2 




C 
































51 


Hnhap3 






T 
















A 




C 




T 


T 




90 


Hnhap4 






















A 




C 




T 


T 




72 


HnhapS 


T 








C 


C 


C 




C 


T 


A 


T 


C 




T 


T 


T 


27 


Hnhap6 


T 






A 




C 


C 


T 






A 




c 


T 


T 


T 




3 



Total 272 



Dot (.) represents the identical nucleotide with the haplotype Hnhapl . 



to calculate a Hd value of 0.662 and a Pi value of 0.0093. 
On the other hand, the yellow Korean native chicken (KY) 
breed had a comparatively low number of variable sites of 
10 and 4 haplotypes, from which a Hd and Pi of 0.591 and 
0.0077 were calculated, respectively. 



A total of 28 nucleotide substitutions were observed 
among the chicken populations which represent 15 
haplotypes (Table 3). Among these, 6 haplotypes were 
identified in the Hyunin population. Of these, two 
haplotypes were unique to the Hyunin population. On the 



Table 3. The identified haplotypes with the mitochondrial D-loop sequence polymorphisms among the Korean native chicken 
populations 



Within 



Nucleotide substitutions 



Total 

sample Haplotype 



popuiauon 
haplotype 


1 


1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


3 


3 


s 


S 


3 


3 


S 


S 


3 


s 


4 






6 


7 


9 


1 


1 


1 


2 


2 


3 


4 


4 


4 


5 


6 


6 


8 


9 


0 


1 


1 


1 


3 


4 


6 


6 


7 


9 


4 


7 


7 


9 


0 


2 


7 


2 


5 


9 


2 


S 


6 


6 


1 


5 


1 


1 


6 


0 


s 


5 


0 


2 


S 


7 


0 


1 


6 


KBhapl 


T 


A 


C 


C 


G 


C 


A 


C 


A 


G 


C 


C 


C 


T 


C 


A 


A 


T 


T 


c 


C 


C 


A 


c 


T 


T 


C 


T 


11 


HI 


KYhapl 


























































7 




Hfhapl 


























































13 




HnhapS 






T 




















































27 


H2 


KYhap3 






T 




















































S 




KRhapl 






T 




















































2 




KBhapS 






T 




















































1 




Hfhap2 






T 
















































A 




10 


H3 


KBhapS 






T 






































T 














2 


H4 


Hfliap4 






T 






































T 














8 




KBhap9 






T 








G 






























T 














1 


H5 


KBhapS 




T 


T 






T 


















T 


G 




C 










G 


T 








C 


4 


H6 


KRhapS 




T 


T 






T 


















T 


G 




C 










G 


T 








C 


S 




KBhap6 






T 






T 








A 












G 
















T 


C 






C 


2 


H7 


KRhap6 






T 






T 








A 












G 
















T 


C 






C 


2 




HfhapS 






T 






T 








A 












G 
















T 


C 






C 


2 




KBhapV 






T 






T 








A 








C 




G 














G 


T 


C 






C 


2 


H8 


KBhapl 


C 




T 






T 




T 






T 




T 


C 










C 


















C 


S 


H9 


KRhapS 


C 




T 






T 




T 






T 




T 


C 










C 


















C 


17 




Hnhap4 


C 




T 






T 




T 






T 




T 


C 










C 


















C 


72 




Hftiap3 


c 




T 






T 




T 






T 




T 


c 










c 


















c 


7 




Yfhapl 


c 




T 






T 




T 






T 




T 


c 










c 


















c 


37 




Hnhap3 


c 




T 


T 




T 




T 






T 




T 


c 










c 


















c 


90 


HIO 


KBhap4 


c 










T 




T 






T 




T 


c 






G 




c 


T 










c 


C 




c 


2 


Hll 


KRhap2 


c 










T 




T 






T 




T 


c 






G 




c 


T 










c 


c 




c 


3 




Hnhap2 


c 










T 




T 






T 




T 


c 






G 




c 


T 










c 


c 




c 


SI 




KRhap4 


c 




T 






T 




T 






T 




T 


c 






G 




c 


T 










c 


c 




c 


3 


H12 


Hnhapl 


c 




T 






T 




T 






T 




T 


c 






G 




c 


T 










c 


c 




c 


29 




KYhap2 


c 




T 






T 




T 


G 




T 




T 


c 










c 


















c 


18 


H13 


KYhap4 


c 




T 






T 




T 


G 




T 




T 


c 




























c 


2 


H14 


Hnhap6 






T 




A 


T 












T 


T 


c 










c 




T 














c 


3 


H15 


Total 




























28 






























437 


IS 



KBhap = Black Korean native chicken haplotype; KRhap = Red Korean native chicken haplotype; KYhap = Yellow Korean native chicken haplotype; 
Hfhap = Hanhyup chicken haplotype; Yfhap = Yellim chicken haplotype; Hnhap = Hyunin chicken haplotype. 
Dot (.) represents the identical nucleotide with the upper haplotype. 
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other hand, four haplotypes were mixed with other Korean 
chicken populations. Of these, four haplotypes (HnhapS, 
Hnhap4, Hnhap2 and Hnhapl) were highly related with the 
KR breed (25 samples) which represents H2, H9, Hll and 
H12, respectively. In the case of the KB breed, only six 
chickens shared three Hyunin haplotypes (Hnhap5, Hnhap4 
and Hnhap2) which are the H2, H9 and Hll, respectively. 
Only one haplotype (Hnhap5) was shared with the KY 
breed (3 samples). Based on the haplotype sharing among 
the Korean native chicken breeds, the Hyunin chicken 
population is very closely related to the KY chicken breed. 
Two chicken populations (Hanhyup and Yellim) had 
Hnhap4 representing the H9 haplotype. On the other hand, 
specific haplotypes for Hyunin population have contained 
33.09% of Hnhap3 (90 samples) and 1.1% of Hnhap6 (3 
samples) which represent HIO and H15 haplotypes, 
respectively. The HIO haplotype, containing the D-loop 
nucleotide substitution C210T, is unique to the Hyunin 
chicken breed. Also, three SNP positions, G212A, C246T 
and C315T, in the D-loop region are specific for H15 
discrimination of the Hyunin breed. Previously, we 
successfully differentiated two commercial chicken 
populations by the selected three D-loop SNPs at positions 
C225T, A239G and C243T (Hoque et al., 2011). Also, it 
was demonstrated that allele-specific SNP typing in mtDNA 
D-loop region at the five SNP sites (C310T, A342G, C446T, 
A686G and C1213T) were SNP haplotypes that can be used 
for identification of the origins of chicken meat (Harumi et 
al., 2011). The D-loop region in mtDNA may also play a 
vital role in the measurement of diversity, and in evaluating 
the population structure of native Korean chickens. 

Network profiling and pliylogenetic analysis 

The estimation of network profiles was traced using 15 



haplotypes in the D-loop nucleotide positions of native 
chicken populations (Figure 1). In our studies, haplotype 
HIO was found to specifically represent the Hyunin 
population that differentiated H9 by the C210T SNP. Also, 
the H15 haplotype for the Hyunin population is 
distinguished from other haplotypes by the three nucleotide 
positions at 212, 246, and 315. Moreover, the H9 Hyunin 
Haplotype is differentiated with H12 by the four nucleotide 
positions at 291, 313, 367 and 370. Besides, the H12 
haplotype is differed from the Hll haplotype by 199 
nucleotide position for the Hyunin haplotype. Only 
haplotype H2 is markedly different from other Hyunin 
haplotypes. These specific SNPs observed in the D-loop 
region will provide valuable markers for the discrimination 
of Korean native chicken breeds. 

The phylogenetic analysis among the chicken 
populations was conducted for identifying relationships 
between Hyunin chicken populations and previously 
identified Korean native chicken breeds (Figure 2). In our 
phylogenetic analysis, haplotypes were clearly 
differentiated among the two groups. Five haplotypes for 
the Hyunin population were clustered together in one clade 
(lower clade in Figure 2) and only one haplotype was 
contained in another clade (upper clade in Figure 2), while, 
three haplotypes in the lower clade were mostly related with 
the KR breed. The internal Hyunin haplotype Hnhap4 is 
also related with two commercial chicken populations and 
the KB haplotype. In addition, Hnhap2 contains the KB 
haplotype. Interestingly, Hnhap3 and Hnhap6 are 
completely separated from other haplotypes, which was 
determined through network profile analysis. On the other 
hand, Hnhap5 in the upper clade is closely related with KR 
and KB haplotypes. Our phylogenetic analyses also gives 
some idea as to the relationship among Korean native 



H_6* 




Figure 1. Network profiles of the 15 haplotypes from the six Korean native chicken populations. The links are labeled by the nucleotide 
positions. The order of the mutations on the branch is arbitrary. The median vectors (mv) are nucleotide junctions. 
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90 



65 

67 



68 



88 



Hfllap4 ] Hanhy up chicken (HI) 



•KBhap9 ^ Blade Korean native chicken (KB) 



gQ K.Bhap5 ] Black Korean natne chicken (KB) 

K.Yhap3 J Yellow Korean native chicken (KY) 

Hfliap2 ]Hanhyupchicl«!n(HO 

I Inliap5 ~\ llvunin chicken (lln) 
KRhap 1 J Red Korean native chicken (KR) 
KKopS 2 BbuMc Korean native chicken (KB) 

KBIiapl ] Black Korean native chicken (KB) 
KYhapl ] Yellow Korean native chicken (KY) 
Hfliapl ] Hanhyup chicken (Hf) 

97|KBhap3 ^ Black Korean native chkken (KB) 
I K Rhap5 Korean native chicken (HQ 



67 



85 



57 



■KBhap7 ^^''"^'''^''■■""■'"''^■'^''''■C^) 



74 



KBhap6 ^ Black Korean native chkiken (KB) 

KRhap6 ~\ Red Korean n;,tiN c chicken (KR) 
HfliapS 3 H^nhyup chicken (lit) 



^— Hnhap6 ] Hvunin chicken fHh) 

43 I KYhap2 

I — KYhap4 



^ cllow Korean nathre chicken (KY) 



67 



84 



49 



KRhap2 ^ Red Korean nathrechxskiai(KR) 
Hnhap2 j Hvunin chicken (Hn) 
KBhap4 3 ^^^^ Korean native chicken (KB) 
K Rhap4 ^ Red Korean native chnksn (KR) 

68'Hnhapl ] nvunit 



96 



70 

66 
66^ 



K miap2 ^ Black Korean native chksken (KB) 

K Kliap3 ^ Red Korean native chksken QSS) 

Hnhap4 ] Hvunin chicken (Hn) 

IIfliap3 ] Hanhyup chwken (Hf) 

Ythapl ]Yellim chksken (YJ) 

HcfaapS jHYtfnintliitlmiftfa) 



0.002 

Figure 2. The neighbor-joining phylogenetic tree using Kimura 2-parameter model with 1000 bootstrap replications in six Korean native 
chicken populations. 



chicken populations. 

In this study, we investigated the D-loop SNPs and 
haplotypes for the identification of Korean native chicken 
populations in a Hyunin private farm. Two previously 
unidentified haplotypes were investigated, and the 
phylogenetic relationships were observed by comparing 
them with haplotypes of other chicken breeds. The results 
indicate that the private farms still hold valuable genetic 
resources which had not been previously identified. To 
maintain this valuable native chicken population, an 
appropriate conservation breeding program is ultimately 
required, and should be conducted with the aid of molecular 
genetic techniques. 
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